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A  THEORETICAL  AND  EXPERIMENTAL  STUDY  OF  STRUCTURES 


IN  THE  MAGNETO-TELLURIC  FIELD 

ABSTRACT 

The  general  problem  of  interpreting  magneto- 
telluric  field  measurements  is  a  difficult  one,  chiefly 
because  exact  solutions  seem  possible  only  for  the 
simplest  geological  geometries.  Because  of  this,  an 
attempt  was  made  to  construct  a  realistic  model  of  the 
earth-ionospheric  current  system.  Measurements  were 
made  on  horizontal  layer  and  fault  models  for  the  case 
of  magnetic  polarization.  These  were  compared  with 
the  exact  solutions  of  Cagniard  and  of  d’Erceville  and 
Kunetz  in  order  to  evaluate  the  model  system.  Agree¬ 
ment  was  found  to  be  reasonable,  considering  the 
complexity  of  conditions  in  the  laboratory. 

The  boundary-value  problem  of  the  dike  in  a  two 
layer  earth  was  then  solved,  again  for  the  case  of 
magnetic  polarization.  Model  measurements  were  extended 
to  include  dikes  and  these  results  when  compared  with 
those  computed  from  the  exact  solution,  were  in  reason¬ 
able  agreement  with  one  another. 
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CHAPTER  I 

INTRODUCTION,  SUMMARY,  AND  HISTORICAL  REVIEW 

The  ever  increasing  interest  in  electro-magnetic 
interactions  with  the  earth  can  be  attested  to  by  the 
number  of  laboratories  pursuing  research  in  this  field 
and  the  number  of  meetings  devoting  time  to  considera¬ 
tions  of  the  subject.  The  subject  is  of  interest  not 
only  from  the  more  basic  scientific  point  of  view  but 
also  from  practical  considerations  in  radio  and  general 
signal  propagation,  and  in  economic  exploration  methods. 

Since  a  major  interest  of  the  geophysics  group  at 
the  University  of  Alberta  is  in  the  field  measurements 
of  the  magneto-telluric  effect,  it  was  felt  that  model 
studies  of  the  phenomenon  would  be  of  considerable 
importance  in  developing  an  understanding  of  the  subject. 
Despite  the  general  importance  of  such  studies,  no  true 
electro-magnetic  model  of  a  layered  earth  has,  to  date, 
been  described  in  the  literature. 

The  difficulties  involved  in  an  electro-magnetic 
model  study  will  become  apparent  in  the  body  of  this  thesis. 
The  difficulties  involved  in  field  measurements  are  already 
well  known  and  well  described  in  the  literature,  and  will 
be  touched  on  at  a  later  point  in  this  work. 
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System  of  Units 

In  all  the  theoretical  developments  to  follow  em 
cgs  units  will  be  employed.  In  all  the  practical  ap¬ 
plications  and  discussion  of  them,  the  so-called  "practical" 
units  will  be  employed.  There  is  no  loss  of  generality 
in  this  procedure  nor  are  there  any  particular  complications; 
this  procedure  merely  follows  that  of  the  important  writers 
on  the  subject  who  are  presently  available  in  English 
translation  and  will  enable  the  interested  reader  of  this 
work  to  read  the  references  with  no  further  transformation 
of  units.  The  "practical"  units  are  those  employed  by 
the  exploration  geophysicist  and  are  related  to  the  em 
units  as  follows: 

H  measured  in  gammas  1  em  unit  =  10^  X 

E  measured  in  mv/km  1  em  unit  =  1  mv/km 

lengths  measured  in  km  1  em  unit  =  10” 5  km 

?  measured  in  ohm-meters  1  em  unit  =  10”H  ohm-meters. 

Thus,  for  example,  skin  depth  which  in  the  rationalized 


mks  system  is 


becomes  in  the  em  cgs  units 


where  T  is  the  period,  and  in  the  "practical"  units: 


Other  important  quantities  have  similarly  simple  trans¬ 
formations  which  will  appear  at  the  appropriate  places 
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in  this  work. 

We  will  consider  only  non-magnetic  media  where  /u  =  i 

and  £  =  eQ  . 

Historical 

All  geophysical  investigations  of  the  subsurface 
nature  of  the  earth1 s  crust  must  of  necessity  depend  on 
contrasts  in  the  physical  parameters  of  the  various 
sub-surface  regions.  One  of  the  most  pronounced  of  such 
contrasts  can  be  observed  in  magnetometer  readings  over 
regions  where  magnetic  materials  lie  included  in  non¬ 
magnetic  country  rock.  Similarly  since  the  resistivity 
of  the  earth  could  be  expected  to  vary  greatly  from 
point  to  point,  one  might  expect  to  observe  pronounced 
contrasts  when  electrical  or  electro-magnetic  measure¬ 
ments  are  carried  out.  An  example  is  the  case  of  a 
saturated  salt  solution  with  a  resistivity  of  the  order 
of  approximately  10"^  ohm-meters  and  dry  salt  with  a 
resistivity  of  the  order  of  100  ohm-meters,  lying  adja¬ 
cent  to  each  other  at  the  edge  of  salt  plugs.  Another 
example  where  pronounced  discontinuities  in  resistivities 
occur  is  the  case  of  massive  sulphides  in  an  unmineralized 
environment.  Indeed  when  the  early  applications  of 
electric  and  electro-magnetic  methods  were  made  they 
appeared  to  offer  tremendous  possibilities  for  successful 
exploration,  and  they  were  and  are  successful  within  the 
limits  which  still  prescribe  them. 
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In  the  earlier  methods  of  exploration  a  source  is 
set  up,  either  D.C.,  in  which  case  the  depth  of  penetra¬ 
tion  depends  on  the  distance  between  electrodes  implanted 
in  the  surface  ol  the  earth;  or  A.C.,  in  which  case  either 
a  magnetic  dipole  or  a  long  line  is  used  and  the  depth  of 
penetration  depends  both  on  the  skin  effect  and  the  ap¬ 
propriate  inverse  geometric  factor.  Whenever  such  sources 
are  used  the  practical  limits  of  exploration  are  restricted 
to  a  few  hundred  feet,  whereas  in  general  the  answers 
required  in  geophysics  have  no  such  modest  requirements. 

Naturally  occuring  telluric  or  ground  currents  were 
first  postulated  by  Davy  in  1821.  Though  Faraday  in  1831 
was  convinced  of  their  existence  as  a  result  of  his  work 
on  electro-magnetic  induction,  the  primitive  equipment 
available  at  the  time  did  not  permit  their  detection. 

Barlow  (1849),  as  a  result  of  studies  of  the  disturbances 
created  in  the  English  telegraph  lines,  was  the  first  to 
obtain  convincing  evidence  that  such  currents  do  flow  in 
the  earth,  these  disturbances  being  most  pronounced 
during  periods  of  magnetic  activity.  Toward  the  end  of 
the  19th  century  many  magnetic  observatories  included 
studies  of  telluric  currents  as  part  of  their  programs. 

The  number  of  stations  doing  so  has  varied,  but  records 
continuous  over  many  years  do  exist.  The  major  preoc¬ 
cupation  of  those  concerned  with  these  earlier  studies 
of  telluric  currents  was  the  origin  and  nature  of  these 


currents 


While  the  correlation  of  telluric  current  and 
magnetic  field  was  appreciated,  the  difficulties  of 
analysis  and  instrumentation  and  the  seeming  lack  of 
positive  results  caused  a  fluctuation  in  interest  in  the 
subject  and  consequently  in  the  number  of  stations  making 
the  dual  observations.  The  diurnal  variations  were 
analyzed  by  Gish  (1939)  on  the  basis  of  only  twelve  world 
wide  stations,  and  on  the  basis  of  his  analysis  he  con¬ 
structed  a  pattern  of  current  distribution  as  shown  in 
Figure  A.  The  increased  activity,  especially  since  the 
inception  of  the  I.G.Y.,  has  resulted  in  a  large  increase 
in  the  number  of  permanent  stations,  forty  as  of  1957, 
with  another  twenty  projected  for  immediate  implementation. 

The  wealth  of  frequency  components  available  makes 
possible  the  study  of  a  wide  range  of  specific  problems, 
ranging  from  frequencies  suitable  for  groundwater,  solid, 
and  liquid  mineral  exploration,  to  geophysical  studies 
of  current  distribution  in  the  lower  mantle  and  core. 

Areas  such  as  Japan  (Rikitake(1952,  1953,  1955)),  and 
Germany  (Bartels  (1954))  display  anomalously  high  magnetic 
variations  with  periods  of  approximately  twenty  minutes. 
This  indicates  anomalously  large  current  distributions 
at  depths  of  some  hundreds  of  km  with  currents  of  the 
order  of  5,000  amperes. 
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One  can  distinguish  between  the  telluric  method, 
which  was  highly  developed  in  France  and  Eastern  Europe 
as  an  exploration  method  shortly  before  World  War  II, 
and  the  magneto-telluric  method,  the  name  coined  by 
Cagniard  ( 1953 ) •  The  possibility  of  the  use  of  the 
telluric  method  for  exploration  was  first  suggested  by 
Schlumberger  in  1921.  In  this  method  a  base  station  is 
used  in  conjunction  with  a  movable  field  station.  It 
obviously  requires  that  the  telluric  currents  are  uniform 
at  leafet  up  to  the  boundary  of  the  region  being  explored. 
We  will  discuss  this  question  in  detail  in  the  next 
section.  This  uniformity  is  all  that  is  required, 
no  further  knowledge  regarding  the  nature  of  the 
phenomenon  being  necessary. 

Some  of  the  workers  in  this  method  considered  the 
associated  magnetic  variations  as  being  the  Biot  and 
Savart  field  due  to  flow  of  current,  in  which  case  E  and 
H  should  be  in  phase,  others  considered  a  rather  simple 
case  of  induction  in  which  case  £  and  H  should  be  9C° 
out  of  phase.  The  scarcity  of  data  on  the  simultaneous 
measurement  and  lack  of  knowledge  of  the  electrical  pro¬ 
perties  of  the  earth  did  not  assist  a  clarification  of 
the  ideas  regarding  the  phenomenon.  However,  as  in  many 
geophysical  methods,  a  lack  of  understanding  of  the 
basic  theory  does  not  preclude  the  possibility  of  some 
success  in  its  application;  all  that  is  really  necessary 
is  uniformity  of  the  source  and  a  contrast  that  can  be 
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measured . 


The  magneto-telluric  method  consists  of  a  simulta¬ 
neous  measurement  of  the  electric  and  magnetic  field 
components  of  an  electro-magnetic  wave;  the  magnetic 
field  measured  at  the  surface  of  the  earth  and  the 
electric  field  measured  as  a  potential  gradient  in 
the  surface  of  the  earth.  The  method  may  be  alternatively 
called  a  "generalized  impedance  method"  since  in  the 
analysis  the  electric  field  and  magnetic  field  appear  as 
a  ratio  „  .  This  concept  of  impedance  is  due  to  Schel- 

n 

kunoff  (1938).  The  magneto-telluric  method,  besides 
having  definite  advantages  in  principle,  especially  in 
that  phase  angle  offers  a  second  diagnostic  criterion, 
also  has  a  practical  advantage  over  the  telluric  method 
in  that  no  base  station  is  required,  and,  further,  the 
uniformity  of  the  field  is  required  to  a  much  more 
restricted  degree.  The  basic  ideas  of  the  magneto-telluri 
method  are  a  recent  development  and  were  motivated  accord¬ 
ing  to  Cagniard  (1956)  by  the  lack  of  correlation  between 
the  predictions  of  the  spherical  harmonic  analysis  of 
Schuster  (1889)  and  Chapman  (1919)  and  the  field  measure¬ 
ments  at  observatories,  as  well  as  inconsistencies  of  the 
telluric  method.  In  1950  Tikhonov  suggested  the  method 
as  a  tool  for  exploration  at  great  depths  while  almost 
simultaneously  Kato  and  Kikuchi  (1950)  showed  that  some 


■ 
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measurements  they  made  on  phase  angles  between  the 
magnetic  and  telluric  fields  could  be  explained  on  the 
basis  of  a  two-layered  earth.  A  series  of  papers  by 
Rikitake  (1952),  Tikhonov  and  Lipskaia  (1952),  and 
Lipskaia  (1953)  were  devoted  to  various  aspects  of  this 
phenomena.  Cagniard  (1953)  published  a  very  comprehen¬ 
sive  paper  on  the  magneto-telluric  method  in  which  he 
developed  graphical  methods  of  considerable  usefulness 
for  interpretative  purposes.  As  Cagniard  envisioned  it, 
the  magneto-telluric  method  was  suited  to  the  large  seal 
exploration  of  sedimentary  basins. 

It  is  one  of  the  facts  of  nature  that  the  earth’s 
crust  is  more  distinguished  by  irregularities  than  by 
regularities,  and  while  layered  sedimentary  basins  do 
exist,  of  equal  interest  are  those  very  irregularities 
which  give  promise  of  providing  economically  valuable 
mineral  resources  or  valuable  clues  to  the  geological 
development  of  the  earth’s  crust. 

Neves  (Ph.D.  Thesis  1957)  and  d’Erceville  and 
Kunetz  (1959)  have  attacked  the  problem  of  a  vertical 
discontinuity,  i.e.  a  fault,  in  a  layered  earth  for 
certain  special  cases.  The  method  of  the  latter  will  be 
discussed  in  detail  in  Appendix  IV. 

The  theoretical  problem  of  a  fault  belongs  to  a 
famous  class  known  as  the  "problems  of  the  finitely 
conducting  wedge".  In  the  more  general  case  these 
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problems  have  not  only  remained  unsolved,  but  have  not 
even  as  yet  been  explicitly  formulated.  In  certain 
special  cases,  which  fortunately  are  geophysically 
significant,  approximate  solutions  can  be  obtained. 

It  is  one  of  the  purposes  of  this  work  to  investigate 
the  adequacy  of  these  solutions.  Another  purpose  is 
to  extend  the  theory  to  certain  cases  of  a  dike,  and  also 
to  study  these  cases  experimentally.  The  theoretical 
development  of  .the  latter  cases  appears  in  Chapter  I V. 
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CHAPTER  II 

DESCRIPTION  OF  THE  MAGNETO-TELLURIC  FIELDS 

General  Properties 

Despite  the  very  low  frequencies  and  the  fact  that 
the  resistivity  ?  is  the  only  parameter  involved,  no 
simpler  formulation  of  the  magneto-telluric  effect  is 
possible  than  is  contained  in  Maxwell* s  equations. 
However,  considerable  simplification  follows  from  the 
experimentally  observed  fact  that  the  sources  of  this 
effect  are  such  that  the  fields  appear  to  be  propagated 
as  plane  waves.  The  sources  are  believed  to  be  vast 
sheets  of  current  in  the  ionosphere. 

Wait  (1954)  has  pointed  out  that  if  the  sources 
were  to  be  considered  as  an  arbitrary  system  of  dipoles 
then  the  fields  as  measured  at  the  surface  of  discontin¬ 
uity  in  the  x,y  horizontal  plane  must  be  related  thus: 


i  +  an 


X  =~  4  Tj  m  a;  -  Eju 


where 


. 


is  the  intrinsic  propagation  constant  of  the  lower  medium, 
and 

c  r 

is  the  intrinsic  impedance  of  the  lower  medium,  O'" 
being  its  conductivity. 

The  above  formulation  takes  into  account  the  curva¬ 
ture  of  the  wave  front  and  the  resulting  phase  angles. 
For^-.l  mho  and  ou  corresponding  to  a  period  of  10  sec. 

—  =  35  km 

a  • 

For  this  case  when  Hx,  Hy  do  not  vary  appreciably  over  a 
region  of  35  km,  the  terms  in  -L  can  be  neglected  and 

E-X  =  V.  H  y 

£y  =  Ux$ 

i.e.,  E  and  H  are  perpendicular  to  each  other.  In  this 
case  the  simplified  treatment  of  Cagniard  (1953),  to  be 
treated  in  some  detail  in  Appendix  II,  is  valid.  For 
longer  periods,  which  are  of  importance  in  the  method. 
Wait  suggests  that  such  uniformity  over  distances  of 
35  km  is  unlikely,  since  the  sources  of  the  magnetic 
variation  are  at  heights  of  only  100  km.  However,  Wait’s 
argument  depends  on  the  assumption  that  the  sources  for 
the  longer  wave  lengths  are  localized,  whereas  it  is 


' 
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precisely  for  the  longer  wave  lengths  that  the  fields 
appear  to  be  uniform,  while  it  is  lor  the  much  shorter 
wave  lengths  that  the  sources  may  be  more  localized. 

In  this  latter  case  the  100  km  heights  of  the  sources 
does  not  impose  such  a  severe  limitation  as  for  the 
longer  wave  lengths  and  sufficient  uniformity  is  still 
apparent  in  actual  measurements.  Cagniard  (1954)  insists 
that  localized  sources  exist  only  in  regions  of  industrial 
installations. 

Schlumberger  and  Kunetz  (1946)  and  Kunetz  (1953, 

1954)  have  confirmed  the  similarity  of  the  magneto- 
telluric  fields  from  regions  as  widespread  as  Madagascar, 
France  and  Venezuela.  The  correlation  of  magneto-telluric 
activity  with  sun  and  auroral  activity  would  seem  to  in¬ 
dicate  the  global  extent  of  the  phenomenon.  In  general, 
the  a  priori  assumption  of  the  magneto-telluric  method, 
that  of  plane  wave  propagation,  would  seem  to  be  justified. 

The  appropriate  equation  for  the  propagation  of 
the  magnetic  field  can  be  written: 

=  +*•«•&-  +  £0  . 

We  will  hereafter  consider  that  field  components  of  both 
»  • 

E  and  H  depend  on  time  as  eJa)^  since  any  arbitrary 
spectrum  of  frequencies  can  be  resolved  in  theory  by 
Fourier  analysis  and  in  practice  by  suitable  filtering. 


In  this  case  the  propagation  equation  is  reduced  to: 


V2W  =  4'tf'jG'oO  -£u>2lf| 


II  - 1 


In  air,  where  o  ,  equation  II-l  becomes: 


while  in  a  medium  such  as  the  earth  where  the  second 

i* 

terra  of  equation  1,  which  is  the  displacement  current 
term,  is  negligible,  the  propagation  equation  reduces  to: 


=  4irj<ru)  w 


Recalling  that  for  both  media  a  =  I  ,  and  f  =  X 

9  °  ca 

we  see  that  while  the  velocity  of  propagation  of  the 
wave  in  air  is  C,  the  velocity  of  propagation  in  the 


earth  is  v  -  fife 


For  a  period  of  one  second  the  velocity 


For  resistivities  as  large  as  10^  ohm-meters,  a  quite 
high  value,  the  velocity  in  the  earth  is  still  only 
approximately  2  x  105  m/sec.  Thus,  by  Snell's  Law, 
a  plane  wave,  even  if  incident  at  nearly  grazing  angle, 
would  be  propagated  into  the  earth  nearly  normal  to  the 
surface.  This  is  a  result  of  the  situation  that  in  a 
conducting  medium  such  as  the  earth,  displacement  current 
is  negligible. 


A  solution  to  the  propagation  equation  for  the 
down-going  wave  in  an  infinite  homogeneous  earth  can 
be  written  in  the  form: 


->  _  4TTj<t”cO 

M  o  c 


^TT<r~oO  z> 


-j  ^TT<5-uJ  z. 


We  define  cl  =  -^==-  in  the  normal  way  as  the  skin  depth, 
the  depth  at  which  the  amplitude  of  the  field  diminishes 
to  —  of  its  value  at  the  surface.  A  list  of  values  of 


e 

the  skin  depth  in  km  versus  periods  and  resistivities 
in  ohm-meters  is  shown  in  Table  I. 

From  Maxwell’s  equation: 


(Vk  H)x  =4'rtvEx 


we  can  write  in  the  case  of  H  polarized  in  the  y  direction 

*?■  M  »• 


-  =  4  IIS"  Ex. 

ZT  *V 


or: 


zT  lEi 
iHv 


In  Table  2  is  shown  the  amplitude  of  H  in  gammas  for 
E=1  mv/km  for  various  values  of  f  and  T  •  The  two  tables 


are  included  at  this  point  under  the  general  description 
of  the  properties  of  the  fields  in  order  to  illustrate 


DEPTHS  OF  PENETRATION  GIVEN  IN  KM  (After  Cagniard) 
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the  possibilities  of  the  magneto-telluric  method,  both 

as  to  the  depth  of  penetration  possible  and  the  easily 

* 

detected  variations  in  the  ratio  of  E  and  H  with  both 
period  and  resistivity. 

The  formulas  of  the  preceeding  paragraphs  are  in 
em  units,  whereas  the  associated  tables  are  in  the 
"practical11  units,  since  for  most  readers  these  are 
easier  to  visualize  both  dimensionally  and  in  magnitude 
than  are  the  em  units. 

In  the  theoretical  development  to  follow,  we 
consider  that  the  surface  of  the  earth  can  be  considered 
as  a  flat,  in  the  dimensions  we  are  considering.  The 
uniformity  of  the  fields  and  the  direction  of  propaga¬ 
tion  in  the  earth  ensure  the  validity  of  the  approxima¬ 
tion. 

For  our  study  we  will  adopt  a  co-ordinate  system 
with  z  positive  downward.  Since  we  will  be  concerned 
with  plane  waves,  which  propagate  normal  to  the  earth’s 
surface  within  the  medium,  we  are  free  to  choose  our 
direction  of  propagation  in  the  positive  z-direction  for 
all  cases.  We  are  further  free  to  orientate  our  x,y 
directions  on  the  surface  of  the  earth  to  suit  the 
direction  of  polarization  of  one  of  the  electro-magnetic 
vectors.  In  practice  E  and  H  will  be  measured  in  mutually 
perpendicular  directions. 
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Let  us  choose  our  co-ordinates  such  that  in  free 
space  at  large  distances  from  the  earth  H=Hy  and  the 
surface  of  the  earth  is  the  x,y-plane  at  z=0.  A 
second  horizontal  interface  when  it  is  introduced  will 
be  at  z=h.  Vertical  discontinuities,  when  they  are 
introduced,  will  lie  in  the  y,z  plane,  i.e.  magnetical¬ 
ly  polarized  cases  only  will  be  discussed  from  a  theor¬ 
etical  point  of  view.  The  surfaces  of  discontinuity 
separate  regions  differing  only  in  their  resistivities. 
Figure  1  shows  the  general  geometrical  arrangements. 

For  simple  layered  systems  with  no  vertical  discontinuities 
we  will  consider  the  surfaces  to  be  infinitely  extended 
in  all  horizontal  directions.  Any  vertical  discontinuities 
will  be  in  the  y,z  plane  and  of  infinite  extent  in  the 
y  direction.  These  may  be  either  dikes  or  faults. 


x 


z-o 


FOFMA  TIOM  / 


Z  =  h 


fz  F OF  M A  T ’/ON  Z 


Fig.  1.  Geometrical  arrangement  of  a  layered  earth 
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Thus  if  we  write  for  two  of  Maxwell’s  equations: 


V  x  M  =  4Ti  + 


it 


V  X  z  = 


-yb 

it 


within  the  media  where  is  negligible  in  comparison 

to  4  sT i  we  can  re-write: 


V  x  M  =  4-TTi  IE -Z  V  x  El  =  -y cub.  IL~3 


From  considerations  of  symmetry  we  see  that: 

( 1 )  All  partial  derivatives  with  respect  to  y 
are  identically  zero* 


( 2 )  Ey  is  identically  zero. 

For  condition  (2)  we  could  as  well  have  written  Bx 
is  identically  zero. 

From  the  second  of  Maxwell’s  equations,  as  written 
above,  we  expand: 


Ak  _^£v 

3y  2)  x 

^  E-Z.  _  ^  Ely 

k  ^1 


b 


£ 


and  observe  that  Bx  and  Bz  and  hence  also  Hx  and  Hz 
are  everywhere  zero  by  the  two  conditions  given  previously. 
We  originally  chose  our  y  axis  such  that  H=Hy  for  z  o 
and  now  we  see -that  W  =  Hy . 
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From  the  first  of  Maxwell's  equations  we  re-write 
and  expand: 


4  UY 


3 


^  W, 


J 


-  41^ Lx 

=  4T^Lz 


which  gives: 


L 


X 


-? 

4T 


r  _  ?  sthly 

4  57  ^  x 


11-4 


IT  -s 


as  a  result  of  condition  1  above. 

From  equation  II-3  we  can  readily  deduce  that  the 
tangential  component  of  E  is  continuous  across  any  sur¬ 
face,  which  giVes  in  our  geometry  that  Ez  is  continuous 
across  a  fault  surface. 

We  observe  at  this  point  that  the  above  deductions 
are  valid  for  the  case  where  one  medium  is  infinitely 
conducting,  ?  =?  o  . 

Applying  equation  II-2,  at  the  upper  surface,  sepa¬ 
rating  regions  where  ?  is  finite  or  zero  from  air  where 
?  is  virtually  infinite,  we  have,  iz  being  zero,  that: 

4Jiy  =  o 

^  X  * 

i.e.,  at  the  upper  surface  Hy  is  constant  even  across  the 
trace  of  the  horizontal  discontinuity  in  ?  . 
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If  we  now  solve  the  propagation  equation  for  It  for 
the  media  we  have  considered: 

V2  W  =  4-li<Taji4 

we  can  then  deduce  the  value  of  it  from  equations  4  and  5» 

Special  Properties  Applicable  to  Model  Studies 

In  Figure  2  are  illustrated  the  relationships  between 
the  various  field  components  across  a  fault  separating  two 
media  in  one  of  which  the  resistivity  is  zero.  In  the 
medium  to  the  left  the  incoming  wave  propagates  downward 
and  may  be  partly  reflected  at  each  interface.  In  the 
medium  to  the  right  where  the  incoming  wave  is  totally 
reflected,  Ex  for  the  incoming  wave  just  cancels 
Ex  for  the  reflected  wave  to  give  a  resulting  value  EX=G. 


r ig .  2 • 


The  magneto-telluric  field  components 
at  a  fault. 
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The  resulting  Hy  which  is  continuous  across  the  strike 
need  not  be  zero  on  the  surface,  since  Hy  for  the  reflected 
wave  is  in  the  same  direction  as  Hy  for  the  incoming  wave. 
Across  the  fault  Ez  is  continuous  and  zero;  i.e.,  there  is 
only  an  x-component  of  E  and  it  is  non-zero  only  on  the 
side  where  P  is  not  zero.  Charge  distributes  itself  on  the 
interface  in  such  a  manner  that  Ex  and  Hy  are  unchanged  from 
the  values  appropriate  to  an  infinitely  extended  case  in 
region  one  and  also  so  that  Ex  is  zero  in  region  two, 
except  possibly  for  edge  effects  in  the  immediate 
neighbourhood  of  the  upper  edge.  If  at  some  finite  distance 
to  the  left  we  introduce  another  discontinuity  separating 
a  third  region  of  resistivity  zero,  extending  infinitely 
to  the  left,  the  same  conditions  will  hold,  therefore 
in  the  region  in  the  centre  which  has  finite  conductivity 
the  field  components  are  those  appropriate  to  an  infinitely- 
extended  layered  system.  In  our  model  we  have  attempted 
to  take  advantage  of  these  properties  and  this  behaviour 
of  the  field  components  in  order  to  simulate  an  infinitely 
extended  media  within  a  rather  restricted  space. 

In  conducting  media  of  finite  resistivity  we  will 
assume  that  in  general  Hy  is  a  function  of  x  and  z; 

E,  since  it  is  derivable  from  Hy,  is  also  a  function  of 
x,  z.  We  have  already  observed  that,  at  distances  remote 
from  vertical  discontinuities,  the  dependence  on  x  must 


vanish  and  we  have  further  observed  how  a  perfect 
conductor  causes  the  field  components  to  be  modified. 

It  is  only  in  proximity  to  structure  in  media  of  finite1’ 
resistivity  then,  that  we  expect  to  find  some  dependence 
of  Hy  on  x  and  therefore  a  z  component  of  E  through  the 
equation: 


£ 


L 


=  JL  4My 
*+>T  x.  • 


EZ  vanishes  on  the  upper  surface  where  Hy  is  a  constant 
due  to  the  fact  that  no  current  flows  across  this  surface 
At  the  vertical  interface  E  must  make  an  abrupt  change 
in  direction  since  Ez  is  continuous  whereas  Ex  is  dis¬ 
continuous.  The  magneto-telluric  ratios,  for  this  case, 
measurable  along  the  upper  surface  where  Ez“0,  is  the 
subject  of  a  paper  by  d’Erceville  and  Kunetz  which  is 
summarized  in  Appendix  IV. 

We  recapitulate  the  qualitative  results  at  this 
point:  at  large  distances  from  the  fault  E  is  in  the 
x  direction,  near  the  fault  it  has  a  z  component,  and 
in  the  intervening  regions  it  has  intermediate  directions 
If  we  now  introduce  a  thin  layer  of  a  perfect 
conductor  on  the  interface,  the  direction  of  E  must  be 
changed  at  least  at  the  surfaces  of  this  conductor  in 
order  that  there  be  only  an  x  component  of  E.  The  lines 
of  force  are  distorted  from  the  previous  case  as 


»• 


illustrated  in  Figure  3-  If  the  conducting  layer  is 
sufficiently  thin  in  comparison  to  distances  over  which 
individual  measurements  are  made,  we  will  assume  that 

9 

the  effect  of  this  infinitely  conducting  layer  can  be 
considered  as  negligible. 


Fig.  3.  Distortion  of  a  field  due  to  a  thin 
layer  of  a  perfect  conductor. 


The  application  of  these  latter  results  to  the 
construction  of  model  faults  and  dikes  will  be  discussed 


later  in  this  work. 


25 


CHAPTER  III 

GENERAL  METHOD  OF  SOLUTION  AND 
INTERPRETATION  FOR  MULTI-LAYERED  EARTH 


Plane  Waves  Incident  on  a  Multi-Layered  Earth 

With  the  co-ordinate  system  previously  described 
and  illustrated  in  Figure  4 


?Z.  FOlZMATlOhJ  Z 


?  h  p&.zma  77  cw  n 


Fig.  4.  Geometrical  arrangement  for  an 
n-layered  earth. 

we  consider  a  multi-layered  conducting  earth  infinitely 
extended  in  the  x  and  y  direction,  the  lowest  layer  of 
resistivity  V>i  extending  infinitely  downward. 

The  equation  for  the  propagation  of  H  in  media  where 
A-  I,  =  degenerates  into: 


= 


-t  X.  -* 
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as  a  result  of  the  form  of  the  time  dependence  and  the 
neglect  of  displacement  current  compared  to  conduction 
current.  For  these  layers  infinitely  extended,  we  can 
further  reduce  this  equation,  by  considerations  o f  sym¬ 
metry,  to  the  ordinary  differential  equation: 

-SLi^y  =  4Tri<ruj^  .  W  -  6 

A  solution  which  satisfies  the  boundary  conditions  for 
any  value  of  x  and  y  and  vanishes  for  z  >co  can  be 
written  down  immediately. 

W*l  =  AL  4-  P'L  ITT  -  7 

for  all  layers  i  «  1,  2,  .  n. 

p-  _  I  _  - 

is  7^  times  di,  the  skin  depth.  The  terms  with  Ai 

Vx. 

represent  an  upgoing  or  reflected  wave  and  the  terms  with 
Bj_,  a  downgoing  wave.  Hence: 

Mn  ^  &n  <2-1 7 i  m  -5 

since  An  must  be  zero  in  the  infinitely  downward-extending 
lowest  layer. 

From  the  continuity  of  Hy  and  Ex  at  each  of  n 
boundaries  we  can  obtain  2n  equations  from  which  we  can 
solve  for  2n-l  constants  A^,  in  terms  of  H0,  the  value 
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of  Hy  at  the  upper  boundary.  This  leaves  one  equation 
with  which  to  solve  for  E0,  the  value  of  Ex  at  the  upper 
boundary.  The  solutions  for  certain  special  cases  have 
been  developed  by  Cagniard.  The  development  is  summarized 
in  the  Appendices. 

Single-Laver  Earth  Infinitely  Extending  Downward 

We  quote  the  following  results  from  Appendix  I: 


and  observe  that  in  this  simple  case  the  electric  field 
leads  the  magnetic  by  45°  and  that  the  resistivity  of 
the  layer  can  be  simply  derived  from  the  magneto-telluric 
ratio . 

Double-Layered  Earth: 

(The  upper  layer  of  resistivity^,  to  depth  h  , 
the  lower  layer  of  resistivity  ^  infinitely  extended 
downward . ) 

The  results  as  derived  in  Appendix  2  give: 


rr  vi  ejc£+<t>-4n 


V2T  N 


where  M,N,<|i  ,  and  ^  defined  in  the  appendix.  Recalling 
from  the  previous  section  that  for  a  single  homogeneous 
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layer: 


£x 

My 


we  observe  that  the  simplest  evidence  of  a  hetereogeneous 
earth  is  the  departure  of  the  phase  angle  from  2+5°. 

In  order  to  use  the  graphical  technique  so  beautifully 
developed  by  Cagniard  (1953)  we  follow  his  procedure  and 
define  the  apparent  resistivity  "  M  through  the  equation: 


Then  returning  to  the  first  equation  of  this  section, 
we  can  write: 


If  one  plots  versus  T  for  each  value  of  f  ,  then  one 


would  obtain  a  separate  family  of  curves  for  each  value 
of  and  h. 

The  phase  angle  which  is  also  diagnostic  of  layer 
ing  is  given  by: 


0  -  ±  -t-  4>  -  y  . 


The  details  can  be  obtained  from  the  appendices. 
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Since  there  is  a  vast  saving  in  labour  to  be  achieved 
by  a  suitable  graphical  method,  we  will  digress  at  this 
point  to  discuss  the  law  of  similitude  and  return  to  des¬ 
cribe  the  method  of  presenting  the  totality  of  curves 
for  the  two-layer  case  in  one  master  set,  immediately 
afterward . 

The  Law  of  Similitude  (Stratton  ( 194 1 ) ) 

Consider  two  geometrically  similar  structures  with  a 
dimensional  ratio: 


L 


^ ,  1 


For  our  electro-magnetic  similitude  we  will  consider: 


^  —  /U  ~  \ 

i 


£  —  £  ==■  to 


9'  =  kff 


The  only  other  variable 


T'  =  ktT 


Similitude  requires  that: 


V  1—1  —  u>  l~l 


implies 


From  the  former  equation 


we  obtain  the  latter,  if: 


V-t  =  lCj>  . 


Further,  from  Maxwell’s  equation: 


V  *  E 


6t  * 


We  see  that  dimensionally: 

C.  =  L 

W  T  ’ 

then: 

_  k  p  E”. 

w  ~  m  w" ,  . 

and 

-  ,A  LS  _ 

^  2T(E)»  ‘‘k?- 

Thus  the  ratio  of  the  apparent  resistivities  is  the  same 
as  the  ratio  of  the  actual  resistivities. 

The  measurements  made  on  models,  as  far  as  size, 
resistivity  and  frequency  are  concerned,  have  immediate 
application  to  a  full-size  earth.  One  need  only  take 
into  account  the  dimensional  ratios 

We  will  proceed  now  to  the  presentation  of  our 
theoretical  curves,  as  for  example,  the  double-layered 
earth  of  the  previous  section. 
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Presentation  of  Data 


The  values  computed  from  the  formula 


can  be  plotted  as  the  logarithmic  ordinates  versus  j It 
as  the  logarithmic  abcissa  for  the  various  values  of 
the  parameters  f{  ,  and  h.  If  one  chooses  ft  =  1  ohm 


meter,  h=l  km,  then  there  is  only  one  parameter  and 
one  obtains  a ’family  of  curves,  one  for  each  value  of  . 


The  complete  set  of  magneto-telluric  curves  is  now 
represented  by  this  one  family,  since  we  can  consider 
the  values  of  f  and  h  as  being  that  of  a  model,  and 
thus  the  curves  obtained  by  plotting: 


can  be  brought  into  co-incidence  with  one  of  the  family 
of  theoretical  curves  by  a  translation  parallel  to  the 
ordinates  of  magnitude  Kp  and  parallel  to  the  abcissa 
of  magnitude  y  .  Referring  once  again  to  Appendix 

2,  one  may  write  alternatively: 


where: 


with  d-[  the  skin  depth  in  the  upper  medium.  The  family 
of  curves  for  has  an  infinite  number  of  points  common 
to  each  curve.  These  are  given  by: 


or 


Z  h 
cL\ 


for  all  n. 


The  point  furthest  to  the  right  has  an  abcissa: 


which  gives  the  result: 

?a  =  ?, 


fr-  8 


for  the  ordinates  of  the  point. 

It  is  customary  to  plot  the  phase  angle  0  as 
the  ordinate  on  a  linear  scale  versus  j~[  as  abcissa  on 
a  logarithmic  scle. 

We  will  finally  transform  to  the  ’’practical"  units 
for  plotting  all  the  results.  In  this  system 


and  the  common  point  of  the  family  of  curves  is  given  by 


One  may  further  observe  that  for  T  sufficiently 
large,  each  curve  of  the  family  becomes  asymptotic  to  a 


' 


finite  value  except  in  the  case  where  =  O  or  oo  . 

This  asymptotic  value,  of  course,  corresponds  to  the 
situation  where  the  upper  layer  becomes  negligible,  i.  e. 
the  skin  depth  is  very  much  larger  than  h.  Similarly 
for  high  frequencies,  the  curve  approaches  the  value  fj  , 
but  only  by  an  infinite  number  of  small  amplitude  oscil¬ 
lations  about  Jj*. 

The  amplitude  of  the  shift  parallel  to  the  abcissa 
that  must  be  simultaneously  made  in  this  translation 
process  in  order  to  bring  the  experimental  curve  into  co¬ 
incidence  with  the  theoretical  curve  is: 


but  by  definition  Kt  is  the  ratio  of  the  abcissa  of  the 
cross-over  point  of  the  theoretical  curve  as  seen  on 
the  experimental  curve: 


T' 

w 


M 


h  =  l 


* 


thus  the  true  thickness  of  the  upper  layer 

10  T~ 

& 

where  ky  is  the  translation  parallel  to  the  axis  of 
ordinates  and  Tl  is  the  abcissa  of  the  cross-over  point. 

y,  —  x  f 

9^  is  the  asymptote  of  the  experimental  curve  for 

r 

large  T  as  read  off  the  extrapolation  of  the  experimental 
curves.  The  interpretation  is  satisfactory  only  if  the 


same  translation  parallel  to  the  abcissa  is  required  to 
bring  the  phase  angle  curve  into  coincidence  with  its 
theoretical  curve.  These  families  of  curves  are  illus¬ 
trated  in  Figures  5  and  6. 

Interpretation  in  Case  of  Three  Formations 

Under  certain  favourable  conditions  the  magneto- 
telluric  method  is  actually  much  simpler  of  interpreta¬ 
tion  than  the  conventional  resistivity  method. 

Consider  a  layered  system  of  three  formations  of 
resistivity  ?,  ,  and  ,  and  of  thickness  h1  ancj  h 

If  7V  is  sufficiently  large  then  it  is  possible  that  the 
h 

effect  of  the  substratum  starts  to  be  appreciable  only 
after  the  effect  of  the  upper  layer  has  ceased  to  be  ef¬ 
fective.  In  other  words,  for  periods  sufficiently  long, 
?a. .  has  approached  closely  to  ?z  and  9  has  approached  4^ 
before  the  effect  of  the  third  layer  has  started  to  make 
itself  felt.  The  curves  of  Figure  7  are  computed  from 
the  theory  in  Cagniard  (1953).  We  can  see  that  if  the 
curves  are  cut  where  indicated,  each  part  can  be  fitted 
to  one  or  another  of  the  family  for  the  two-layer  case. 

Such  favourable  conditions  do  not  always  exist, 
but  a  good  deal  of  interpretation  can  be  carried  out  with 
a  minimum  catalogue  of  computed  curves. 


Apparent  Resistivity  in  Ohm-Meters 
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Fig.  5.  Master  curves  of  apparent  resistivity  for 
magneto-telluric  soundings  over  a  two- 
layer  earth. 
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Fig.  6.  Master  curves  of  phase  angle  for  magneto- 

telluric  soundings.  (after  Cagniard,  1953) 


©  ( 


Computed  curves  for  three-layer  case. 

h*  =  900  h,  f,  ,Pt  and  f3  in  the  ratio  2:10:1. 

(after  Cagniard,  1953) 


CHAPTER  IV 


EXTENSION  TO  THE  THEORY  FOR  DIKES  IN  A  LAYERED  EARTH 


In  Appendix  IV  are  derived  the  magneto-telluric 
ratios  for  a  fault  according  to  dfErceville  and  Kunetz 
(1959)*  We  have  summarized  there,  only  the  case  where 
the  vertical  discontinuity  extends  from  the  upper 
surface  through  a  single  layer  to  the  top  of  a  uniform 
infinitely  downward  extending  lower  layerf 


tz 


Fig.  8.  A  fault  in  a  double-layered  earth. 

The  fault  surface  is  the  y-z  plane  at  x*=Q  and 
extends  infinitely  in  the  positive  and  negative  y 
directions.  (Figure  8). 


Of  at  least  equal  interest  to  the  fault,  from  a 
geological-geophysical  point  of  view,  is  the  case  of  a 
dike.  For  this  reason  it  was  felt  advisable  to  add  to 


the  already  existing  theory,  solutions  for  some  interest* 
ing  cases  of  this  type. 

Consider  a  dike  consisting  of  two  vertical  planes 
infinitely  extended  in  the  ±  y  direction  separating  a 
region  of  resistivity  ,  from  the  medium  on  either  side 
of  resistivity  . 


Z  =  h 


Y; 


Fig.  9.  A  dike  in  a  double-layered  earth. 


We  take  our  origin  such  that  the  yz  surfaces  of  dis¬ 
continuity  lie  at  x  3  ,  4  •  Media  (1)  and  (2)  are  both 

underlain  by  a  common  infinitely  downward  extending 
substratum.  (Figure  9)« 


At  large  distances  from  the  discontinuity  the 
fields  would  be  those  of  the  appropriate  undisturbed 
double  layer.  The  solutions  for  these  cases  are  worked 
out  in  the  Appendix  III.  The  solution  for  H  in  the  un¬ 
disturbed  region  in  the  medium  of  resistivity  /,  will 
be  designated  as  Hj  and  that  in  the  medium  with 
as  H2,  bearing  in  mind  that  H=Hy.  We  will  now  assume 
that  the  total  solution  for  the  field  H  is  a  function 
of  x  and  z  of  the  form: 


in  region  1,  with  a  similar  expression  with  subscript  2 
in  region  2.  Each  term  of  this  sum  must  satisfy  the 
propagation  equation  independently. 

We  will  assume  that  we  can  write  the  disturbance 
term  as  the  product  of  two  functions: 


We  will  further  assume  that  g(z)  can  be  expanded  in 
a  Fourier  series. 

We  will  consider  two  special  cases: 

Case  I:  Infinitely  Resistive  Substratum 

In  Appendix  III  appear  the  solutions  for  a  simple 
double  layer  of  this  type.  We  can  observe  that  the 
solution  H^(z)  vanishes  for  z=h.  Reference  to  this 
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appendix  will  indicate  that  Hy (x, z )=Hi ( z )  +  Pi(x,z) 
must  also  vanish  for  z=h.  Furthermore  we  have  seen  in 
Chapter  II  that  Hyis  continuous  for  z=0  even  across  the 
trace  of  discontinuities  in/.  Since  P(x,z)  is  assumed 
to  have  no  effect  for  large  x  we  must  have  that 
P(x,o)=G.  Thus,  since  P(x,z)  is  a  function  which  vanishes 
for  z=0  and  z=h,  we  can  expand  the  term  g(z)  into  a  sin 
series  of  argument  in  the  interval  o,h  and  include 

the  coefficient  of  sin  in  the  f(x)  term. 

Pl  fx,z)  =  1  fm  (x)  M,r)  r\T(x  * 


Each  term  of  this  series  must  satisfy  the  same  propaga 
tion  equation: 


<3ZU  , 

this  implies: 


>-+-  S  C  y  c~  £0  i— I 

'  t 


n(x)  ^un  n'STz 
h 


m2J(  2‘ 


-f  4  T TjCTcd 


fin  Of? 


n  Sc  z 


which,  after  cancelling  the  common  sin  term  has  the 
solution: 

fin  C  _k|n 4  bin  ek'** 

where: 


) 


' 


:  •  ' 
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This  solution  satisfies  the  condition  that  it  vanish  for 
|X|  cx=>  if  we  take  only  the  a-type  terms  for  positive 
x  and  only  the  b-type  terms  for  negative  x  both  in 
medium  2,  while  both  positive  and  negative  exponentials 
are  permissable  in  the  finite  region  of  medium  1. 
Symmetry  conditions  permit  us  to  equate: 

=  k>t.n 

which  we  observe  gives  the  same  value  for  Hy(x,z)  for 
the  same  positive  and  negative  values  of  x.  Utilizing 
the  symmetry  condition  is  entirely  equivalent  to  using 
one  of  the  boundaries  for  solving  for  continuity  of  the 
field  components,  only  the  remaining  boundary  is  now 
available  for  solving  for  the  remainder  of  the  coef¬ 
ficients  and  we  write: 


-  ^ 1  n  ^ 

&  ~T]r 


n=i 


I  " 


fcinl 

zTT” 


—  CL 


n  u  z 

h 


2. 


If  we  now  expand  H2-  H-^  into  a  sin  series  of  the  same 
argument: 


^2.-^1  ~  ^  ^  n 

n  =  ! 


n  7T  z 

h 


we  can  equate,  at  the  boundary,  term  by  term  to  obtain: 


i  n 


^jiA. 

2  h 


a 


\  Z  h  / 


W 


Recall: 


r  =S-  ^A.y 
4»r  J7 


and  equating  term  by  term  at  x  =  L 
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a*ne 


~^z.n  £ 
Z  h 


O  TV  -  i  o 


Solving  for  a(n  from  equations  (<=? )  and  (10) : 


I  M 


A 


_  fj  n _ 

Z  4-  km  *»£  'UTlAfEir^I 

r27T  kznhA  Izh 


E- 


a 


■2  K7 


-  An  i  4^  w.  i^ftr) 

K 


^  hii  co^  +  4<,n6 


4tn 


*  h 


£  h 


IST-  iz 


Following  the  well  known  procedure  of  Fourier  series 
expansion,  the  coefficient  of  sin  irf( z  in  the  expansion 
of  H2  **  H1: 


A„  -  ij(U4-U 


where : 


,-M,  =A,<2^  +  B.e1^  *-Ae^  -B.e  ^ 


Z  'I 


Pi 


L  /4-?T^ 


uJ 
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quoting  the  result  of  Appendix  III,  Equations  A-3  and  A-4 


A;. = 


HL) 


g>.  —  -  Up'*- 


>{J  ht 


2  i  hL 


from  which: 


A.  = 


ZVj  n  ^  ^  "  hf  ■ -) 

'  k  2  k 


i  r? 


substituting  back  into  equations  (11)  and  &)  we  obtain: 


CL\n  = 


2  "tr;  fh,z-h,2W  Mo  n. 


t  2 1-  2 

*•  i  n 


2  ('fini')  +  hikL  2  U  fJiini  ^ 
L  V  2h  y  iThT  k  jzh  /J 


^2n 


2Tj  (h/-  k,2)W0  n«  oink  (%£) _ 


k  2  k 

K  i  n  ** 


i  n 


In  medium  I 


Uv/  = 


u, 


y  h,)|_ 


_/SG,ni+'(i%+<s}n'h'  ftp, 


ex? 


+)air12a«fi^  ^ 


n:  | 


Ex= 


U„S' 


O  J] 


4-7 T  Z  'x*dr\ (y  j  h,) 


v7  , 1  0  * ' 

Lr. e  +  r  -  J 


Cc 

?,  W 


4R"  j_  H  a'n 


cOMi  CCl  11(1 


h 


h 


n=  I 


At  the  upper  surface  z  =  0,  where  we  measure  the  magneto- 
telluric  fields  Hy  =  H0 
and : 


i  n  i 


where: 


U 


4 

L 

n»l 


h2  (li|Kj  h  ) 


k|n  k 


Zn 


cwH(  V)+ 

v  2k  /  k^hz  Uh  / 


kzMh,: 


In  medium  2 


The  constant  term  is  the  same  as  for  medium  1  with 
the  subscripts  interchanged  but  the  disturbance  term  is: 

cy? 


K 


V  -k2M*  n7Tz 

-  \  dZn  <2  ~7T 


n»l 


and  the  disturbance  term  for  Ex: 
no 


h  =  I 

At  the  surface  where  Hy  =  H0  : 


U 


z 


where : 


0^5 


L 

n-  i 


k2Mt  ( £  \ 

£  ~iyr  $vdr\  (  2h  j  e 


-fcjLn* 
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k,„n 


k 
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Case  II:  Infinitely  Conducting  Substratum 

In  Appendix  III  apear  the  solutions  for  a  simple 
double  layer  of  this  type.  The  solution  ( z )  does  not 
vanish  for  z  =  h.  Thus  in  the  general  expression: 

U(x,z)  =  (z)  +  Pl  (*,z) 

P^(x,z)  must  vanish  for  z  =  0  for  the  same  reason  as 
in  the  preceeding  example,  while  at  z  =  h  it  must  have 
a  finite  value  determined  by  the  continuity  of  H(x,z) 
at  .  Thus,  neither  nor  P  can  be  rep¬ 

resented  simply  by  a  sin  series  in  the  interval  0  to 
h.  For  convenience  we  will  develop  these  functions  in 
a  sin  series  of  argument  ,  where  in  order  to 

satisfy  the  boundary  conditions  at  z  =  h  while  vanishing 
at  0  and  2h  we  will  define: 

Hi- U,  “  f(2h-z),  btZ^2.b. 

Proceeding  as  in  the  previous  example: 

P,  4W1  yp  ^ 

each  term  of  which  satisfies: 

-  4  (TjV  jj  (4 

results  in: 


' 


which  has  solutions: 


r  ~  ^  I  r7  y  b  1  h.  V 

tin  =  a,'n  e  Th  +  am  e  ^  t 
where  in  this  case: 

km  =  |  n2^2,  +  4jln{Z  t 

and  we  have  equated  a^n  =  b^n  lor  the  same  reasons  as 
in  the  previous  example.  The  equations  of  continuity 
for  Hy  and  Ez  at  x  =  -{-  result  in  the  same  solutions  lor 
ain  in  terms  of  An  except  that  h  is  replaced  by  2  h: 


noting  that  in  Kqn  and  K2n,  h  is  also  replaced  by  2  h. 
When  solving  for  An  we  note  that: 


A 


r  ,  h 


f(z)  $\n 


V\  m 


j.  n 


Jli  + 


h  -z)bi 


b\ n 


yifTz 

Zh 


di 


J 


where : 

f  (z)  =  H-W, 

and  H2  and  are  obtained  in  Appendix  III. 

A„=  i<=  Kj  JJ0(VA2)  . 

K,n  Kzn 
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i.e.,  only  the  odd  terms  of  the  series  appear,  therefore 
we  must  substitute  2n+l  for  n  in  K^n  and  K.2n»  with  the 
understanding  that  n  starts  at  zero  instead  of  at  unity. 

The  complete  result  for  an  infinitely  conducting 
substratum  is: 


U, 


cosh  [  k,  2h  ) 


k,n  k  2 


'in 


similarly : 


Co 


5^1 


fh 


+ 


Kin 


4  h 
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Presentation  of  Data  and  Interpretation  ior  Faults  and 

Dikes 

The  records  of  a  magneto-telluric  traverse  will 
in  general  contain  a  broad  band  of  frequencies  recorded 
at  each  station.  While  a  frequency  analysis  and  inter¬ 
pretation  may  be  attempted  at  an  individual  station, 
using  the  methods  outlined  in  Chapter  III,  such  a 


. 


procedure  could  lead  to  erroneous  conclusions  if  it  were 
based  on  an  assumption  of  a  layered  sedimentary  basin. 

It  would  be  more  meaningful  to  construct  and  analyze 
profiles  at  discrete  frequencies.  These  profiles  would 
be  obtained  by  plotting  the  amplitudes  and  phase 
angle  respectively  for  the  appropriate  station  co¬ 
ordinates.  Because  of  the  large  contrast  in  resistivity, 
and  hence  magneto-telluric  response,  it  is  useful  to 

p 

plot  the  amplitudes  of  g  as  a  logarithmic  ordinate. 
However,  for  obvious  reasons,  positions,  regardless  of 
the  units  chosen,  must  be  plotted  on  a  linear  scale. 

The  theoretical  expressions  for  the  magneto-telluri 
ratio  become  increasingly  complex  as  the  number  of 
interfaces  is  increased,  even  for  a  simple  layered  system 
When  there  are  vertical  discontinuities  it  is  easily  seen 
that  it  is  no  longer  possible  to  represent  the  totality 
of  magneto-telluric  soundings  by  a  single  family  of 
curves,  as  in  the  case  of  the  one  interface  of  the  two¬ 
layered  system. 

The  interpretation  for  the  case  of  a  fault  as 
discussed  by  d’Erceville  and  Kunetz  (1959)  illustrates 
the  order  of  the  increase  in  the  number  of  families  of 
curves  required  for  a  complete  interpretation.  The 
addition  of  a  second  vertical  interface  for  the  dike 
adds  another  parameter,  namely  nJin  the  width  of  the  dike, 
and  hence  another  order  of  families  of  curves. 
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We  will  restrict  our  efforts,  and  those  of  the 

University  of  Alberta  computing  centre,  for  the  most 

part  to  the  construction  ol  those  curves  from  which  we 

can  attempt  to  interpret  our  experimental  data. 

It  can  readily  be  seen  that  the  first  term  in  the 

theoretical  expressions  for  ^  is  just  the  result  one 

h 

obtains  for  the  undisturbed  double-layer  cases  of 

Appendix  III.  The  summation  for  region  2  which  lies 

outside  the  dike,  converges  very  rapidly  for  £>0  and 

h 

indeed  is  virtually  zero  lor  x=h.  Thus  one  can  use 
the  method  of  interpretation  of  Chapter  III  to  solve 
for  Jz  ,  and  h  from  the  spectrum  ol  frequencies,  at 
stations  sufficiently  far  removed  from  the  disturbing 
boundary.  In  this  case  fz  refers  to  the  resistivity 
of  the  upper  layer.  Within  the  dike  the  range  of 
may  be  very  restricted  and  hence  one  is  in  general 
unable  to  solve  for  as  in  the  case  of  fx  . 

In  the  summation  term  the  multiplying  factor, 
besides  constants,  contains  parameters  in  the  following 
combination: 


for  the  region  inside  the  dike,  and: 


for  the  region  outside  the  dike.  We  have  plotted 
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£  versus  the  ratio  £  for  the  various  values  of  the  ratio 
H  h 

*£  and  of  h  and  £  which  have  been  used  in  the  model. 

7 

We  have  transformed  our  data  into  the  "practical"  system 
of  units  for  this  purpose.  The  theoretical  and  experi¬ 
mental  results  will  be  described  and  discussed  in  the 
next  chapter. 

The  theoretical  results  obtained  by  d»Erceville 
and  Kunetz  (1959)  for  the  fault,  indicate  certain 
features  of  interest  which  may  be  applicable  also  to 
the  dike.  Some  computed  curves  for  this  case  will  be 
shown  in  Figures  12  and  13,  Chapter  VII.  Besides  the 
expected  result  that  at  long  distances  from  the  dis¬ 
continuity  has  the  value  appropriate  to  the  undisturbed 
layered  system,  is  the  interesting  result  that  the  ratios 
of  the  values  of  ^  on  adjacent  sides  of  the  discontinuity 
are  just  the  ra tics  of  the  D.C.  values  of  f  . 

The  aforementioned  appropriate  value  of  §  at 

H 

remote  distances  depends  on  two  factors  in  addition  to 
’ the  resistivities  of  the  upper  layer.  The  resistivity 
of  the  substratum,  and  the  period,  can  dominate  the 
behaviour  of  for  all  positions  measured. 

The  curves  of  Figures  12  and  13  illustrate  the 
remarkable  effect  of  the  substratum  for  longer  periods. 

It  can  be  seen  that  the  curves  labeled  A  are  flat  right 
up  to  the  discontinuity  on  both  sides,  whereas  curves 


. 


B  display  large  variations  near  the  interface.  For 
the  former  case  it  can  be  seen  that  only  at  the  highest 
frequency  used  and  then  only  in  the  region  of  lowest 
resistivity  of  the  upper  layer,  is  there  a  beginning  of 
a  departure  from  a  constant  value. 
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CHAPTER  V 

THE  REAL  EARTH  AND  THE  LABORATORY  MODEL 

One  of  the  difficulties  involved  in  the  measurement 
of  natural  electro-magnetic  fields  is  the  extremely  com¬ 
plex  nature  of  the  earth’s  electrical  properties.  Some 
values  for  resistivity  are  shown  in  Table  3  from  Rooney 
(1935): 


Media  Ohm-meters 

Clays,  dense  alluvia,  etc .  1  to  30 

Sedimentary  rocks,  new  .  1C  to  300 

Sedimentary  rocks,  old  .  300  to  2,000 

Igneous  rocks  .  500  to  5,000 

Coarse  gravel,  sand,  etc .  750  to  50,000 


Table  3  Resistivity  of  Various  Earth  Constituents 

According  to  the  method  of  interpretation  discussed 
previously,  it  should  be  possible  to  ignore  relatively 
thin  layers  of  highly  conducting  materials  and  cor¬ 
respondingly  thicker  layers  of  less  conducting  materials 
in  the  event  that  the  appropriate  range  of  frequencies 
is  chosen.  In  the  earth  an  upper  layer  of  sedimentary 
rock  which  may  be  as  thick  as  10  km  is  underlain  by  a 


, 


53 


granitic  layer  to  a  total  depth  of  30  to  2+0  km.'  It  can  be 
seen  from  the  table  that  the  resistivities  of  these  two 
layers  overlap,  hence  they  could  be  considered  as  a 
single  layer  underlain  by  a  basaltic  layer  of  resistivity  10^. 
Computations  as  reported  by  McNish  (1939)  give  a  value  of 
10  ohm-meters  to  a  depth  of  100  km  with  indications  that 
below  that  the  resistivity  is  1  ohm-meter.  It  is  obvious 
that  these  values  do  not  compare  with  any  of  the  known 
values  as  given  above  and  must  then  be  the  effect  of  a 
highly  conducting  region  corresponding  to  a  depth  pene¬ 
trated  by  the  twenty-four  hour  diurnal  variation.  The 
vertical  discontinuities  can  be  more  pronounced  and 
closer  to  the  surface.  A  survey  of  the  available  data 
illustrates  the  severe  restrictions  on  the  possible 
materials  suitable  for  the  construction  of  a  scale  model, 
bearing  in  mind  that  the  scale  factors  must  satisfy  the 
relationship, 

K l  =  Kj>  Kt  . 

Apart  from  solutions  the  only  homogeneous  conducting 
materials  are  graphite  and  the  metals,  both  of  which  would 
require  inconveniently  small  models,  and  hence  impossibly 
severe  restrictions  on  the  size  of  probe  spacing  and 
magnetic  detectors.  From  among  the  solutions  we  chose 
that  of  ordinary  NaCl,  since  it  has  a  reasonably  low 
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resistivity  in  concentrated  form  and  it  has  a  relatively 
slow  reaction  rate  with  certain  useful  metals.  The 
temperature  coefficient  of  resistivity  is  also  quite 
reasonable  for  use  in  an  open  room.  The  resistivity  of 
the  concentrated  salt  solution  at  22°  C.  is  4  x  10“^  ohm- 
meters.  The  salt  solution  provides  a  uniform  single 
layer  which  is  in  our  laboratory,  underlain  by  a  small 
depth  of  a  dielectric  and  then  by  concrete  floor  under¬ 
neath  which  are  metallic  conduit  and  pipe.  The  details 
of  the  model  are  given  in  the  next  section.  For  a  dike 
we  have  used  graphite  which  gives,  even  with  the  concen¬ 
trated  solution,  an  extreme  contrast  from  a  geological 
point  of  view.  We  have  taken  advantage  of  the  assumed 
behavious  of  thin  layers  of  very  highly  conducting 
material  to  separate  solutions  of  different  concentration 
and  hence  conductivity  as  was  discussed  in  Chapter  II. 

This  has  allowed  us  to  simulate  dikes  and  faults 
of  moderate  resistivity  contrast.  The  details  of  these 
model  structures  will  be  discussed  in  the  next  chapter. 

Another  major  problem  of  field  measurements,  that 
of  suitably  recording  low  amplitude,  low  frequency 
signals  against  a  persistent  background  of  higher  frequency 
diurnal  variations,  is  being  attacked  from  an  instrumental 
point  of  view,  as  is  the  equally  formidable  problem  of 
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resolving  the  various  frequency  components,  on  the 
records,  for  analysis.  A  laboratory  model  avoids  this 
complication,  since  in  the  experiment  one  has  control 
not  only  of  the  set  of  discrete  frequencies  one  wishes 
to  use,  but  also  of  the  geometry  and  thus  the  direction 
of  polarization.  The  possibility  of  anisotropic  resist¬ 
ivity  variations,  which  have  been  discussed  by  Cantwell 
(I960),  can  add  considerably  to  the  difficulties  of  both 
field  and  laboratory  studies. 

While  a  model  study  can  greatly  simplify  certain 
aspects  of  the  problem,  especially  in  verifying  theor¬ 
etical  solutions,  there  are  special  problems  in  electro¬ 
magnetic  modeling,  in  addition  to  those  which  plague  all 
other  model  studies.  These  problems  arise  due  to  the 
coupling  of  remote  objects  with  the  model  through  the 
electro-magnetic  field.  In  this  way  end  effects  are 
far  more  severe  than  in  the  electro-static  case,  for 
example.  Radio  stations  and  beacons  can  be  quite  readily 
detected.  Methods  of  reducing  this  latter  effect  without 
recourse  to  filtering  will  be  described  in  the  next  chapter. 
The  effect  of  unwanted  conductors,  some  of  which  are  also 
magnetic,  is  a  problem  of  the  laboratory  in  which  the 
experiment  is  carried  out.  In  our  laboratory  the  conduit 
and  service  pipes  under  the  floor  provided  the  most 
serious  problem  in  that  it  is  no  longer  permissible  to 
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speak  of  a  simple  two-layered  problem  in  the  discussion 
of  our  model. 

In  order  to  appreciate  the  effect  of  the  finite 
dimensions  of  the  tank,  consider  a  conducting  medium 
threaded  by  an  alternating  magnetic  field.  The  induced 
electric  field  will  cause  a  current  to  flow,  and  if  the 
medium  is  finite  in  length  in  the  direction  of  the  current, 
charge  will  pile  up  in  such  a  way  as  to  partly  neutralize 
the  field. 

Preliminary  measurements  in  our  model  indicated 
that  the  magnitude  of  the  magneto-telluric  ratios  were 
many  orders  of  magnitude  too  small  to  verify  the  theory 
for  a  layered  earth.  In  addition,  the  magnitude  of  E  was 
so  small  that  the  probe  spacing  was  of  necessity  made 
too  large  to  bring  out  the  details  of  the  effect  of 
vertical  dikes  or  faults.  We  were  in  fact  measuring 
average  fields  over  distances  which  covered  the  whole 
range  of  variations. 

Since  our  major  purpose  was  to  investigate  the 
effect  of  vertical  discontinuities  we  adopted  certain 
modifications  to  our  model  in  an  endeavour  to  allow  the 
use  of  smaller  probe  spacing  and  to  diminish  the  end 
effects.  These  modifications  are  discussed  in  detail  in 
the  next  chapter. 


CHAPTER  VI 


APPARATUS  AND  PROCEDURE 

The  magnetic  field  was  generated  by  a  current 
which  was  as  large  as  60  amperes,  flowing  in  a  rect¬ 
angular  loop  constructed  of  £"  x  l/6M  copper  bus  bar. 
(Figure  10).  The  power  was  supplied  by  a  McIntosh  60- 
watt  audio-amplifier  matched  to  the  loop  by  a  step-down 
transformer  with  a  single  turn  secondary.  The  audio¬ 
oscillator,  power  amplifier  and  frequency  counter  are 
shown  in  Figure  B.  The  turns  ratio  could  be  adjusted 
from  4:1  to  40:1  in  discrete  steps.  The  loop,  with 
dimension  3m  x  6m,  lay  in  a  vertical  plane  with  the 
nearest  long  side  50cm  above  the  centre  line  of  a  flat 
wooden  tank.  (Figure  C) .  This  centre  line  was  the 
x  axis  of  our  co-ordinate  systems  and  was  parallel  to 
the  current  flow  in  the  long  sides  of  the  loop.  The 
tank  contained  a  salt  solution  which  acted  as  the  con¬ 
ducting  upper  layer  of  ouf  model.  Dikes  and  faults  lay 
with  their  strikes  perpendicular  to  the  x  axis, 
i.e.  in  the  direction  of  the  magnetic  field.  The  tank 
had  a  length  of  2m  and  a  width  of  lj  m. 

Preliminary  measurements  indicated  that  the  end 

effects  were  extremely  pronounced.  The  E  ratio  was 

H 

many  orders  of  magnitude  too  small,  while  at  the  same 
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tig.  10.  Jchematic  diagram  ol  general  experimental  arrangement. 
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Fig.  B.  Instrument  section;  audio-oscillator, 

power  amplifier,  and  frequency  counter. 


Fig.  C.  Wooden  tank,  containing  layer  of 
solution  and  graphite  dike. 
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Fig.  D.  Photographs  of  signals  on  oscilloscope  screen. 
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time  the  ratio  varied  by  over  7$  in  a  traverse  of  1  50  cm 
from  the  centre  of  the  tank  and  this  effect  increased 
radically  as  the  ends  were  approached  more  closely. 

We  have  attempted  to  improve  our  model  by  placing 
copper  sheets  at  each  end  of  the  tank  and  shorting  them 
together  with  the  bus  bar.  This  bus  bar  is  led  out  at 
right  angles  to  the  current  flow  and  joined  2h  meters 
away  at  the  wall  of  the  laboratory  (Figure  10),  in  order 
to  minimize  the  interaction  of  the  current  in  this  part 
of  the  circuit  with  the  current  in  the  conducting 
solution.  The  resulting  E.M.F.  induced  in  this  loop 
did  not  appear  to  affect  the  measurements  too  seriously. 
The  uniformity  of  the  field  along  the  x  axis  was  greatly 
improved  while  the  ratio  appeared  to  be  of  the  right 
order  of  magnitude. 

The  measurement  of  the  magnetic  field  was  ac¬ 
complished  by  using  a  pick-up  coil,  previously  calibrated 
in  terms  of  the  free  space  response  of  the  coil  with 
respect  to  the  computed  magnetic  field  at  the  centre  of 
a  circular  loop. 

The  electric  field  was  derived  from  the  voltage 
drop  between  two  probes  separated  by  a  distance  of  2  cm 
along  the  x  axis.  It  is  thus  the  average  field  which  was 
measured,  over  a  distance  which  was  made  as  short  as 
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feasible.  In  actual  practice  a  third  probe  was  intro¬ 
duced  which  was  connected  to  the  centre  tap  of  the 
measuring  instrument,  in  our  case,  a  Techtronix  502 
oscilloscope.  In  this  way  our  measurement  of  E  is: 

E  .  -  Va  fob  -  VflMVq-Vq) 

b-a,  b-a 

where  Vg  is  the  voltage  to  ground. 

Both  the  magnetic  and  electric  signals  appeared 
in  turn  as  the  vertical  component  of  a  Lissajous  figure 
on  the  screen  of  the  oscilloscope,  the  horizontal  com¬ 
ponent  being  a  reference  signal  which  was  measured  as 
the  voltage  drop  across  a  non-inductive  element  of 
the  source  loop.  The  measured  and  computed  resistance 
of  this  element  was  U»&U  x  10”^  ohms. 

The  third  probe  was  introduced  into  the  electric 
field  measuring  system  in  order  to  cancel  out  that  part 
of  the  noise  which  affected  both  probes  alike.  This 
function  it  performed  most  efficiently;  the  noise,  which 
appeared  when  an  attempt  was  made  to  measure  simply 
Va-Vg,  was  many  times  larger  in  magnitude  than  the  signal 
level,  but  could  be  reduced  to  tolerable  levels  when  the 
differential  probe  system  was  used.  That  part  of  the 
noise  which  was  directional  and  had  a  component  in  the 
x  direction  could  not  be  so  eliminated.  Our  conducting 
layer  was  apparently  an  efficient  "aerial”  for  radio 
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frequencies  which  were  radiated  upward  through  the  floor. 
This  high  frequency  noise  which  modulated  all  electric 
signals  could  be  "stopped"  in  a  visual  presentation  on 
the  oscilloscope.  The  tuning  requirement  was  tbo  exact¬ 
ing  to  allow  a  photograph  to  be  taken  of  the  modulated 
signal  and  in  general  this  noise  merely  produced  a  line 
broadening.  It  was  felt,  however,  that  the  signal  to 
noise  ratio  was  adequate  without  the  use  of  filters. 

The  magnetic  pickup  was  not  as  susceptible  to  noise  since 
its  response  was  peaked  in  the  audio  range  of  frequencies. 

Both  the  electric  and  magnetic  signals  were  photo¬ 
graphed  on  Polaroid  film.  The  lens  gave  a  1:1  image  of 
the  oscilloscope  screen,  the  oscilloscope  itself  being 
previously  calibrated.  Triple  exposures  were  made  to 
facilitate  the  measurement  of  the  horizontal  and 
vertical  amplitudes  and  the  intercepts,  each  individual¬ 
ly.  Sample  photographs  appear  in  Figure  D. 

In  order  to  simulate  a  fault,  a  tray,  one  end  of 
which  is  a  thin  copper  plate,  is  placed  in  the  tank, 
the  copper  plate  then  separates  a  region  containing  a 
concentrated  salt  solution  of  resistivity  4  x  10  ohm- 
meters  from  a  region  containing  a  dilute  salt  solution 
of  resistivity  .71  ohm-meters.  The  depth  of  the  layer 
was  20  cm  in  this  case. 
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In  the  case  of  the  dike,  a  tray  of  length  10  cm 
with  copper  plates  on  both  ends  was  used.  Two  situations 
were  examined;  the  one  in  which  the  dike  contained  the 
concentrated,  more  highly  conducting  solution,  in  an 
environment  of  the  less  conducting  solution;  and  the 
situation  with  the  solutions  interchanged . 

While  it  is  assumed  that  it  is  permissable  to 
neglect  the  effect  of  thin  metallic  conductors  where  the 
current  and  electric  field  lines  are  normal  to  the  metal¬ 
lic  surface,  it  is  not  permissable  to  introduce  such 
conducting  material  which  provides  conducting  paths 
parallel  to  this  field  direction.  Hence,  in  order  to 
study  discontinuities  lying  at  depth,  it  is  necessary 
to  use  some  other  device,  than  surrounding  a  solution 
with  a  metallic  conductor.  Two  solid  dikes  were  used, 
extending  10  cm  in  the  x  direction  as  in  the  previous 
cases  but  of  half  the  depth.  This  would  correspond  to 
a  large  depth  of  burial  in  the  real  earth,  depending  on 
the  frequency  employed.  A  graphite  dike  provided  a 
sharp  low  resistivity  contrast  and  a  block  of  concrete 
encased  in  wood  provided  a  sharp  high  resistivity 
contrast,  in  each  case  compared  to  both  a  concentrated 
and  dilute  salt  solution. 

For  comparison,  traverses  were  run  on  these  same 
dikes  in  a  shallow  layer  of  the  concentrated  solution. 


In  this  case,  2  mm  of  the  solution  overlay  the  dikes  in 
order  for  the  electric  probes  to  make  contact. 

The  results  and  interpretation  for  all  these  cases, 
as  well  as  various  magneto-telluric  soundings  on  the 
Cagniard  system  will  be  discussed  in  Chapter  VII. 


CHAPTER  VII 


RESULTS,  INTERPRETATION,  AND  CONCLUSIONS 


In  our  calculations  for  the  theoretical  curves, 
we  have  employed  the  actual  dimensions,  frequencies 
and  resistivities  which  were  used  in  the  experiment. 

To  these  ’’real”  measurements  correspond  an  infinite 
number  of  ’’models”  of  earth  size  dimensions,  according 
to  the  scaling  factors  we  choose  for  any  pair  of  the 
parameters;  the  third  parameter  being  determined 
through  the  relationship: 

=  kf  Kt. 

In  Table  4  are  shown  several  possible  geological 
situations  corresponding  to  the  several  examples 
selected  from  among  our  measurements. 

Frequency  Analysis  of  the  Layered  Earth 

In  the  two-layered  earth  case  it  can  be  observed 
that  at  the  shortest  period  used,  the  upper  layer 
thickness  was  less  than  one  skin  depth.  According  to 
the  method  of  interpretation  discussed  in  Chapter  III 
for  a  simple  layered  earth,  it  is  shown  that  for  suf¬ 
ficiently  short  periods,  the  apparent  resistivity 
approaches  the  resistivity  of  the  upper  layer,  or 
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phrasing  it  more  simply,  the  skin  depth  is  such  that 
the  fields  virtually  do  not  penetrate  the  lower  layer. 


In  Figure  11  are  shown  the  results  of  plotting 
from  the  measurements  made  on  the  two-layer  model: 

iP-  -  ll  I 


to  the  same  scale  as  the  master  curves  of  Figure 
These  curves  are  extrapolated  to  shorter  periods  to 
indicate  their  possible  asymptotic  approach  to  the 
value  of  resistivity  measured  in  a  DC  conductivity 
cell . 


At  the  longer  period  end  it  can  be  seen  that  in 
the  three  cases  presented,  the  apparent  resistivities 
coincide.  For  these  very  long  periods  we  are  essential 
ly  measuring  the  common  resistivity  of  a  lower,  relativ 
ly  higher  conductivity  medium.  We  could  estimate 


f  <  I05  . 

We  have  superimposed  on  these  experimental  curves 
segments  of  the  master  plot  for  the  two-layer  case. 

It  will  be  recalled  from  the  discussion  on  interpreta¬ 
tion  in  Chapter  III  how  a  three-layer  curve  can  be 
accurately  derived  in  certain  cases  as  a  composite  of 
two-layer  segments.  Under  less  restrictive  conditions 
a  composite  curve  can  indicate  qualitative  aspects  of  a 
more  general  layered  earth. 
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•  —  h  =  i  x  io 
o  —  h  =  2  x  io-4 
a  —  h  =  2xio-4 


Fig.  11.  Frequency  analysis  for  a  simple  conducting 
layer  measured  in  the  tank. 


We  interpret  the  results  plotted  in  Figure  11  as 
indicating  that  the  relatively  thin  upper  layer  of  our 
model  is  underlain  by  an  extremely  thin  layer  of  more 
highly  resistive  material,  which  is  underlain,  in  turn, 
by  a  relatively  thick  layer  of  highly  conducting  material. 
Keeping  in  mind  that  relative  thickness  is  in  terms  of 
skin  depth,  one  can  see  that  a  few  millimeters  of  a 
conducting  metal  can  be  much  thicker  for  our  purposes 
than  a  few  meters  of  concrete.  Our  model  is  in 
reality  a  complicated  multi-layered  system.  The  salt 
solution  is  in  a  wooden  tank,  for  all  practical  purposes 
an  extremely  thin  layer  of  dielectric,  which  in  turn 
rests  on  an  uneven  concrete  floor,  thus  there  are  also 
air  gaps.  Imbedded  in  the  concrete  is  the  usual  struc¬ 
tural  steel  and  beneath  that,  more  concrete,  service 
pipes,  conduit,  and  earth.  A  quantitative  analysis 
of  this  arrangement  was  not  considered. 

Discussion  of  Theoretical  Results  for  the  Dike 

In  Figures  12  -to-  17,  are  plotted  the  theoretical 
curves  for  faults  and  dikes  with  the  parameters  as  listed 
in  Table  5.  The  experimental  results  for  these  cases  are 
also  plotted  on  the  same  figures.  The  computations  which 
have  been  completed  by  the  computing  centre  to  date, 
have  parameters  selected  to  fit  the  experimental  data. 

Even  at  the  highest  frequency  used,  the  periods  are 


Figure  Curve  Type  of 

No.  No.  Structure 
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PARAMETERS  OF  COMPUTED  CURVES  PLOTTED  IN  FIGURES  12  to  17. 
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Fig.  12.  Computed  and  experimental  values  of  magneto- 
telluric  ratio  over  a  fault  at  the  higher 
frequency . 
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P, 


Fig.  13.  Computed  and  experimental  values  of  magneto- 
telluric  ratio  for  a  fault  at  the  lower 
frequency . 
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Fig.  14.  Computed  and  experimental  values  of  the  magneto- 
telluric  ratio  for  a  conducting  dike  at  the 
higher  frequency. 
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Fig.  15.  Computed  and  experimental  values  of  the 
magneto-telluric  ratio  for  a  conducting 
dike  at  the  lower  frequency. 
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Fig.  16. 


Computed  and  experimental  values  of  the 
magneto- telluric  ratio  for  a  resistive 
dike  at  the  higher  frequency. 
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Fig.  17*  Computed  and  experimental  values  of  the 

magneto-telluric  ratio  for  a  resistive  dike 
at  the  lower  frequency. 
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inconveniently  long  for  a  complete  interpretation. 
Computation  is  continuing  at  the  present  time  in  order 
to  remedy  this  situation. 

Even  so,  certain  marked  features  can  be  observed. 

In  direct  contrast  to  the  case  of  the  fault,  it  is  pos¬ 
sible  to  have  the  ratio  of  |j  on  the  adjacent  sides  of  a 
dike  interface  differ  from  the  DC  case.  For  long  periods 
two  possibilities  exist.  For  an  infinitely  resistive 
substratum,  the  situation  is  dominated  by  the  substratum 

and  the  discontinuity  in  H  is  the  same  as  that  of  the 

H 

DC  case  at  each  interface.  For  an  infinitely  conducting 

substratum,  the  discontinuities  themselves  seem  to 

dominate  the  situation  and  £  at  one  interface  is  affected 

H 

to  some  extent  by  the  more  distant  interface. 

For  shorter  periods  the  situation  is  more  complex. 

In  addition  to  the  effect  of  either,  or  both,  of  the 
horizontal  and  vertical  interfaces,  there  is  the  effect  in 
the  upper  layers  of  the  skin  depth. 

The  curves  of  Figures  12  through  17  indicate  that, 
for  the  frequencies  we  have  used  to  date,  the  upper  medium 
reacts  like  a  linear  resistor,  the  current  being  constrained 
to  flow  in  relatively  thin  upper  formation  by  the  zero 
conductivity  of  the  lower  formation. 

The  same  condition  is  not  true  for  the  zero 


resistivity  substratum.  Currents  are,  in  this  case,  free 
to  flow  unrestricted  by  geometrical  boundaries.  Hence, 
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even  though  the  depth  of  the  upper  region  is  negligible, 

E 

to  the  extent  that  at  remote  distances  tt  would  be  the 

H 

same  in  both  regions,  pronounced  effects  in  the  region 
about  the  discontinuity  can  be  observed. 

The  striking  similarity  between  the  correspond¬ 
ing  curves  of  Figures  12  through  17,  which  differ  only 
in  the  frequencies  used,  is  probably  only  useful  as  a 
criterion  that  both  frequencies  are  relatively  low. 

The  same  information  is  obtainable  in  the  case  of  a  fault 

E 

from  the  fact  the  values  of  are  the  same  at  remote 
distances  from  the  fault  in  either  media. 

The  phase  angle  variations,  for  the  frequencies 
which  we  have  measured  and  computed,  are  not  sufficiently 
pronounced  to  warrant  treatment.  Further  results  from 
the  computing  centre  will  be  more  useful  in  the  study 
of  phase  angles,  since  higher  frequencies  will  be 
employed . 


Experimental  Results 

The  measured  values  of  £  are  plotted,  together  with 

H 

the  appropriate  theoretical  results,  in  Figures  12  through 
17.  It  can  be  seen  that  the  magnitudes  of  the  experimen¬ 
tal  results  are  in  reasonable  agreement  with  curves 
corresponding  to  a  zero  resistivity  substratum.  This 
is  also  seemingly  in  accord  with  our  previous  conclu¬ 
sion,  based  on  the  frequency  analysis,  that  our  complicated 
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model  of  a  layered  earth  was  underlain  by  a  very  low 
resistivity  substratum.  It  should  be  pointed  out  that 
an  extrapolation  of  concepts  based  on  a  two-layered  case 
to  a  more  complicated  case  is  hazardous.  We  have 
already  seen  that  the  conditions  of  the  horizontal 
layering  can  have  a  very  pronounced  effect  on  the 
response  due  to  a  vertical  discontinuity. 

'We  have  mentioned  in  Chapter  VI  that  we  measure 
the  average  field  over  the  distance  of  the  probe  spacing. 
This  is,  of  course,  the  same  as  the  procedure  in  actual 
field  measurements.  Further,  in  our  model  we  are  unable 
to  measure  right  up  to  the  interface  and  hence  are 
unable  to  measure  in  the  most  sensitive  region  for 
detecting  variations  in  the  field  components.  Despite 
these  facts,  the  results  do  indicate  a  curvature  in  the 
correct  direction. 

In  Figure  18,  we  have  plotted  the  results  of 
measurements  in  which  a  graphite  dike  of  resistivity 
8  x  10“"  ohm-meters,  and  a  resistive  dike  of  resist¬ 
ivity  larger  than  1G°  ohm-meters,  and  whose  lateral 
dimensions  are  10  cm  square,  have  been  immersed  in  a 
concentrated  salt  solution  of  depth  20  cm  and  resistivity 
U  x  10  ohm-meters.  If  we  consider  here  scaling  factors 
of  10^  then  these  correspond  to  dikes  of  depth  1  km  buried 
to  a  depth  of  1  km  of  overburden.  In  the  case  of  the 
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Fig.  18.  The  magneto-telluric  ratio,  on  a  linear  scale, 
over  high  contrast  dikes  with  deep  overburden. 


graphite  dike  this  could  correspond  to  massive  sulphides 
surrounded  by  country  rock,  whereas  the  resistive  dike 
could  be  considered  as  an  igneous  intrusion.  We  have 
plotted  the  magneto-telluric  ratio  on  a  linear  scale 
because  of  low  contrasts  observed  in  the  signals  which 
have  been  attenuated  by  the  1  km  of  conducting  material. 

In  Figure  19,  are  plotted  the  results  obtained 
from  traverses  across  these  same  two  dikes  immersed 
in  a  shallow  layer  of  concentrated  salt  solution  just 
deep  enough  to  provide  sufficient  solution  over  the 
dikes  to  make  electrical  connections  with  the  probes. 

This  layer  of  depth  2  mm  would  correspond  to  20  m  of 
weathering  using  the  same  scaling  factor  as  in  the 

i* 

previous  mentioned  case. 

Though  the  possibility  of  detecting  either 
conducting  or  non-conducting  bodies  at  appreciable  depths 
is  of  considerable  interest,  the  theoretical  development 
for  a  geometry  of  this  type  has  not  yet  been  attempted. 
However,  some  understanding  of  the  effect  of  overburden 
can  be  obtained  from  a  study  of  these  curves. 

In  Figure  20,  is  shown  the  variation  in  the 
magnetic  field  over  a  conducting  dike.  We  have  plotted 
this  on  a  linear  scale,  in  addition  to  the  logarithmic, 
since  the  effect  is  quite  small. 
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Fig.  19*  The  magneto-telluric  ratio  f<?r  high  contrast 
dikes  with  shallow  overburden. 


H(X)  LOG  SCALE 


Fig.  20.  The  variation  in  the  magnetic  field  over  a 
discontinuity . 


H  iy)  LINEAR  SCALE 
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CHAPTER  VIII 

CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  WORK 


The  value  of  a  model  study  lies  primarily  in 
its  role  of  advancing  the  understanding  of  phenomenon 
not  directly  susceptible  to  mathematical  analysis. 

Before  it  is  possible  to  proceed  with  confidence  in 
applying  model  results  to  a  real  earth,  quantitative 
verification  of  cases,  for  which  analytic  solutions  do 
exist,  is  essential.  While  the  experimental  results 
of  this  work  are  in  reasonable  qualitative  agreement 
with  the  existing  theory,  more  exact  correspondence 
would  be  desirable.  To  this  end,  a  model  incorporating 
a  better  approximation  to  either  or  both  of  the  infinite¬ 
ly  conducting  or  infinitely  resistive  substratum  should 
be  developed.  The  use  of  much  higher  frequencies  would 
assist  in  this  project,  since  the  dimensions  of  the 
model  could  then  be  reduced,  if  the  result  of  introducing 
shorted  end  plates  is  as  effective  as  we  assume.  Further 
experimentation  regarding  this  last  point  is  required. 

In  the  cases  discussed  in  this  work,  where  the 
magnetic  field  is  polarized  parallel  to  the  structure, 
the  validity  of  the  assumption  that  the  tangential  com¬ 
ponent  of  the  magnetic  field  is  constant,  is  established 
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by  the  results  shown  in  Figure  20.  We  feel  confident, 
therefore,  in  proceeding  to  complete  a  catalogue  of 
computed  curves.  By  varying  the  different  parameters, 
sets  of  curves  can  be  obtained  which  will  clearly 
illustrate  the  effect  on  the  magneto-telluric  field  of 
various  structures  and  resistivity  contrasts. 

The  solution  to  the  boundary  value  problems 
where  the  electric  vector  is  parallel  to  the  strike  of 
the  structure  would  be  desirable.  In  the  practical 
application  in  exploration  known  as  '•  Afmag" ,  the  varia¬ 
tion  in  the  plane  of  polarization  of  the  magnetic  vector 
is  measured  near  or  at  the  surface  of  the  earth. 

The  vertical  component  of  the  magnetic  vector 
is  believed  to  be  due  to  the  variation  of  the  magnetic 
field  across  discontinuities  when  the  polarization  is 
in  the  direction  as  described  above.  In  a  more  academic 
application,  the  long  period  variations  in  the  vertical 
component  of  the  magnetic  field  as  measured  at  coastal 
observatories  could  perhaps  be  explained  in  terms  of 
polarization  at  depth.  In  any  event,  the  solution  to 
this  type  of  problem  is  inherently  more  difficult  than 
the  cases  described  in  this  work,  since  one  can  no  longer 
assume  the  constancy  of  either  field  vector  in  crossing 
the  trace  of  the  discontinuity. 
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APPENDIX  I 

This  is  a  special  case,  the  geometry  of  which  is 
illustrated  in  Chapter  III,  Figure  U,  where  n  =  1, 
fn-f.  and  the  single  layer,  in  keeping  with  our 
original  description,  extends  infinitely  downward. 

The  propagation  equation: 

-ffiV-  =  4rztncUy 

a  z 

has  a  solution: 

IL  -  Be  ^ 
y 

where  p  has  been  defined  previously  as: 

D  =  _ ! _ 

fPfcrZo  . 

There  is  no  reflected  wave  since  the  medium  is  infinitely 
extended  downward. 

At  the  upper  boundary,  continuity  of  Hy  requires: 

Lip'c)  Mc 

where  H0  is  the  magnetic  field  in  air,  Hy  is  the  medium. 
Therefore:  - 

uY=  *  . 

From  Equation  U,  Chapter  II: 

r  ?  ML 

LX=  *4fr  -jt 
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provides  that: 

4s(  p  ° 

At  the  upper  surface  where  z  = 


0: 


■1°) 


-X  -  4Trp  1  'o  . 

Thus  the  measurement  of: 


fx.  I  ,  j  :  4i(  ~u) 


t-l 


y  io 


\ J  (4Tr)"  <■ 


U  UT 


relates  the  magneto-telluric  ratio  to  the  resistivity 
of  the  medium  and  the  period  T  of  electro-magnetic 
variation.  We  re-write  this  expression: 


/  bx 

R 


y 1  o 


__  iff 

T  J  4 

Uj  e 


Thus,  for  an  infinite  single  layer,  E  leads  H  by  45° 
and  the  resistivity  of  the  medium  can  be  expressed: 

f  =  2Tl^ 


u 


y 


APPENDIX  II 


A  TWO-LAYERED  EARTH 


This  is  a  special  case  of  the  multi-layered  earth 
illustrated  in  Figure  4,  Chapter  III,  with  n  =  2.  Forma' 
tion  2  of  resistivity  fz  extends  infinitely  downward. 

The  solution  to  the  propagation  equation: 


is: 


=  ATT JCTCO  Py 

a.  ^ 


IT. 


V  p7 


( III— 6 ) 


UL  =  AL£  P:  ( III-7 ) 


P  =  — A - 

r  r - 


In  Formation  1  the  solution  is: 

U>  =  A,e  p‘  +5,  e  p| 

while  in  Formation  2  the  solution  is: 


’y 


Ex  is  obtained  from  Hy  by  application  of  equation 
II-4  which  in  Formation  1  is: 


r  =  -A  J-L 
4T r  P, 


4,e  -5,e 


while  in  Formation  2: 
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The  continuity  of  Hy  at  z  *=  h  results  in: 

hz 


A 


e 


vTh,  n  _-(Th2 


+  E>.e  ' 


=  e 


A- 1 


where 


h  -  ~ 

h*  P, 


The  continuity  of  Ex  at  z  =  h  results  in: 


p, 


A  Vj  h  n  "(T  h.l 

■A  £  +  D.e  =  ~  z 


J 


p* 


A-£ 


-i/T  U 

Eliminating  <&  z  between  equation  A-l  and  A-2: 


/TP  .  n  'U~K  ft  i  a  »5"h 


r  -Ch, 


0 


A,  &  £>0 


-2/Jh, 


^  L 

ygj  -  fe_~ 

^  +  A1, 


-A, 


<2 


which  when  substituted  back  in  equation  A-l  gives 

fj  ( h ,  -  hi) 


6>,  = 


2.  ->•  V r  1 


z^cr, 


c? 


and: 


zj< 


fcr.fc.  e^(^  +  h^ 


Thus,  at  the  surface  z  =  0,  we  can  write: 


E 


M 


*  ___  1  _ 

Uy  fz^f  N 


0 


-L($+4,-‘f') 
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where : 
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iv i  oin 
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NJ  co^>h  Y 


N  Sin  y 
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+  1  sinh  Jl 

cos 

\n 

a, 

a, 

dz 

a.j 

d-i 

Pi 

sinh  il 
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Osh  n 

-sin 

h 

a, 

a. 

r_L 

sinh  il 
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cosh-t 
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COS 

h 

[d-i 

oi, 

oi-z 

dL 

a, 

Z L 

coshil  4-  -L 

sinh  -tl 

sin 

Ji 

d 

i  a 

d-2 

d, 

* 


and  where: 

ciL  = 


is  the  skin  depth. 

If  one  defines  the  apparent  resistivity: 


then  it  can  easily  be  shown  that: 


2h 

9 a_  i  4-cos  eh 

~c~  ~  '  L  o  2  h 

->i  m  +  m  ~  2  cos  -j- 


where: 


m 


The  analysis  and  interpretation  is  done  in 


Chapter  III 
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APPENDIX  III 

SPECIAL  CASE  OF  A  TWO-LAYER  EARTH 
WITH  A  FINITELY  RESISTIVE  UPPER  FORMATION 

We  will  treat  in  this  appendix  both  the  case  of 
an  infinitely  resistive  and  an  infinitely  conductive 
substra turn. 

Case  1:  Infinitely  Resistive  Substratum 

t' 

The  geometry  is  the  same  as  in  the  previous 
example  of  the  general  two-layered  earth  illustrated 
in  Figure  4,  Chapter  III,  with  n  =  2. 

The  propagation  equation  in  this  case  is: 


where  the  second  term  on  the  right  is  included  to  take 
into  account  propagation  in  Formation  2  where  because 

9  only  the  displacement  current  term  is  non-zero. 
In  Formation  2  a  solution  has  the  form: 


since  there  is  at  most  a  downgoing  wave.  It  is  obvious 
that  there  is  no  attenuation  in  Formation  2,  which  is 
what  is  to  be  expected  in  a  non-aissipa tive  medium. 

In  Formation  1,  where  the  displacement  current  is 
negligible : 


1-1,= 

F?- 


37V; 
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and  at  the  boundary  z  =  h, 


A 


IT  ST 


-Jr  h 


Pi  +  B  *  B, 


For  any  reasonable  value  of  h,  the  thickness  of  the 
upper  stratum,  and  for  periods  of  interest  for  the 
magneto-telluric  method,  the  exponent  on  the  right 
hand  side  of  the  equation  is  approximately  zero.  Thus, 
writing  in  the  usual  way: 


h,  -  £ 

1  P,  » 


we  obtain: 


/“A  i  £2 


/TK  +  B.e'^'  =  & 


and  from  the  application  of  Maxwells  equation  written 
in  the  form: 

V  =  4-?r<r£  + 

O  t  , 

we  have  for  our  one  dimensional  geometry: 


aw 


ci  2- 


y 


4  ?T<r  4j£ 


6C> 


£ 


x . 


For  the  upper  formation  where  displacement  current  is 
negligible  we  obtain  Equation  II-U: 

t-,  ^  4TT  JZ  , 

whereas  in  Formation  2  where  6^  =  o 

r  _  J  AiL 

~  <3  2. 


% 


At  the  boundary  where  z  =  h: 

A  ,/n,'-e><2-'Th'1=  -Jb: 


-i  iS 

4rr  p,  L 


J 


Eliminating  B2  between  the  two  equations,  we  obtain: 


;  A,  ■  = 

A,  =  E>, 

[y7I  %  +  1] 


-\J  h 


In  the  approximation  already  introduced,  where  ao 
corresponds  to  periods  of  a  few  seconds  and  when  skin 
depths  are  larger  than  a  few  meters,  we  can  write: 

a  n 

in  which  case,  by  substituting  back  into  the  equation 
relating  B2  to  A|  and  B^ ,  we  find  to  the  approxima¬ 
tion  we  have  made  that: 

=  O  , 

At  the  upper  surface,  z  =  0: 


A.+  B  =  M 


o 


therefore : 


B.  =  l~lo  e 


(T 


2  sinh»(jfj  /i()  , 


A  ~  3 


and : 


11 

A  = 

1  2iinh(/Jh,) 


A -4 
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The  rather  surprising  result  that  the  magnetic 
field  vanishes  at  the  surface  of  an  infinitely  resistive 
substratum  may  appear  at  first  glance  to  be  antithetical 
to  the  case  of  an  infinitely  conducting  substratum,  on 
the  surface  of  which  the  electrical  field  must  vanish. 
However,  the  situation  is  rather  similar  to  the  flow  of 
energy  in  a  wave  guide  in  which  there  is  a  mismatched 
impedance  element.  Thus,  that  portion  of  an  incoming 
wave  which  is  not  reflected  at  the  first  suface  is 
trapped  and  dissipated  within  the  conducting  medium. 

Writing  the  complete  expression  for  the  electric 


field : 


we  obtain  for  the  magneto-telluric  ratio  at  z  =  0: 


A-  5 


Case  2:  Infinitely  Conducting;  Substratum 

Proceeding  as  in  the  forementioned  Case  1,  we  can 

write : 


» 
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then  writing  the  equation  of  continuity  for  Hy  at  z=G, 
and  for  Ex  at  z=h: 


A,  +  B,=  U0 


Eliminating  between  these  two  equations: 


we 


^1 

MoAl^h' 

2-  cosln  f|/j  In,') 

A-6 

obtain: 

Br 

Wo  1 

£  cosh  (fj  h() 

A-7 

r  — 

4  JT 

fi 

p, 

-i CTh.  /Jf-  /7h,  -i/7|- 
2  2  1  -  <2  2 

■ 

J 

I 


and  finally: 


a -a 
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APPENDIX  IV 

A  FAULT  IN  A  TWO-LAYERED  EARTH 

We  will  discuss  two  special  cases  where  the  sub¬ 
stratum  resistivity  is  zero  or  infinite. 

Case  1:  An  Infinitely  Resistive  Substratum 

The  geometry  is  illustrated  in  Figure  8,  Chapter 
IV.  Region  1  of  resistivity  ?{  and  region  2  of  resis¬ 
tivity  ?z  are  underlain  by  a  uniform  formation  of 
infinite  resistivity  infinitely  extended  downward. 

Following  the  procedure  of  d'Erceville  and  Kunetz 
(1959)  outlined  in  Chapter  IV,  we  write  the  solution  in 
the  upper  formation: 

(V,y)  =  Hi  iz)  +  PL  (x^)  ' 

For  either  region,  where  ?i(x,z)  is  the  term  which  is 
added  to  the  solution  for  the  simple  layered  system 
because  of  the  disturbance,  Pj_(x,z)  is  a  function  which 
vanishes  forjx'| — and  also  for  z  =  0  and  z  =  h. 

We  therefore  assume  that  we  can  write  Pi(x,z)  as  a 
product: 

R  (x,z)  =  fi(x)  3i(z) 

where  further,  g±(z)  can  be  expanded  as  a  Fourier  series 
of  argument  .  We  proceed  further  and  include  the 

coefficient  of  sin  ~-r—  with  fi(x)  to  write  finally: 


9o 


(x t y )  -  h-  (z;  +  y  fin  (v)  im  ^ 

n 

where  each  term  on  the  left  must  satisfy  the  propaga¬ 
tion  equation: 


As  in  Chapter  IV  this  leads  to  the  expression: 


tin  (x)  =  am  e  k'°  *  +  bin 

ktn  =  irL+jh~z  ,  h;  =  h  f+TT6l,Z 


In  the  region  of  negative  x  only  the  term  with  a 
positive  exponential  can  be  valid,  thus  only  the  bqn 
coefficients  are  non-zero.  Similarly,  only  the  aqn 
coefficients  are  non-zero  in  the  positive  x  region.  At 
the  boundary  between  regions  1  and  2  which  is  the  y-z 
plane  at  x  =  0,  we  require  that  both  Hy  and  Ez  are 
continuous : 


P  (■ °,z )  +  P  P)  =  P  (0,1.)  +  Hz  (£), 


Thus,  if  we  expand  H2  -  Hi  as  a  Fourier  series  of  the 

same  argument  as  in  Pq-  we  can  equate  at  the  boundary 

Q  \  hi 

term  by  term  both  in  Hy  and  ^  .  The  coef¬ 
ficient  An  of  sin  in  the  Fourier  expansion  of 

H  2  -  Hq  is  found  by  the  standard  method: 


Pn  I  j  ~  Pi)  Sl 


n  dz 

h 


The  boundary  conditions  then  yield: 
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®-in  D; in  - 

^/n  n  +  ^n  ^°2h  ~  c - 


The  expression  ibr  H2  -  Hi  for  this  case  is  found  from 
equations  A-3  and  A-4  of  Appendix  III  and  the  integra¬ 
tion  of  An  is  readily  performed: 


ci 


in 


=  2TjU0(V-  h*) 


n 


k  1  L  1 

r|  +  ft/  .kin.  1 

Hf2  kzn  J 

with  a  similar  expression  for  b2n* 
Writing  H(x,z)  explicitly: 


+  2Sr,‘M o(V-htZ)  ^ 


HI  r 

n<£  sin 


M  H  2 


hikin'; 


1 


we  can  find  Ex  and  finally  the  magneto-telluric  ratio 
at  z  =  0  in  both  regions.  In  region  1: 


u  _ 


^  cotm  (fj  -zrryu,  h  u, 


6^7 


u,  = 


Z  ^m|h 

n 


/ 

h-/ 


k,/  ki/  f|  4  hLhin. 


The  magneto-telluric  ratio  for  region  2  is  found 
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from  the  preceding  expression  by  a  simple  interchange 
of  subscripts. 

Case  2:  An  Infinitely  Conducting  Substratum 

In  this  case  all  the  results  of  the  previous 
section  hold,  except  that  exactly  as  in  Chapter  IV  for 

this  substratum  we  expand  g ( z )  as  a  sin  series  of 

nTTz. 

argument  —  since  in: 

R.  <*>z)  =  fu(*)  9i.fc) 

we  must  have  a  disturbance  function  which  vanishes  for 
|xi-^oo  and  z  =  0,  but  which  has  a  finite  value  at 
z  =  h.  Therefore,  when  we  also  expand  H2  -  Hi,  as 
a  Fourier  sin  series  also  in  we  must  define: 

Mi.- 1-1 , =  F(z)  otz  z  1 

M lz-M,=  E(zh-z)  h~z£zh 

where  H2  and  Hi  are  the  appropriate  solutions  for  the 
simple  double-layer  system  of  Case  2,  Appendix  III. 

When  this  procedure  is  carried  out,  we  obtain 
the  magneto-telluric  ratios: 
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with: 


kin  =  /  (zni  \)ZU'2'+  4j  h,z 


is  again  obtained  from  —  by  a  simple  inter¬ 
ne  Ho 

change  of  subscripts. 
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